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Abstract. We analyse the off-resonant Raman interaction of a single broadband 
photon, copropagating with a classical 'control' pulse, with an atomic ensemble. It 
is shown that the classical electrodynamical structure of the interaction guarantees 
canonical evolution of the quantum mechanical field operators. This allows 
the interaction to be decomposed as a beamsplitter transformation between 
optical and material excitations on a mode-by-mode basis. A single, dominant 
modcfunction describes the dynamics for arbitrary control pulse shapes. Complete 
transfer of the quantum state of the incident photon to a collective dark state 
within the ensemble can be achieved by shaping the control pulse so as to 
match the dominant mode to the temporal mode of the photon. Readout of 
the material excitation, back to the optical field, is considered in the context of 
the symmetry connecting the input and output modes. Finally, we show that the 
transverse spatial structure of the interaction is characterised by the same mode 
decomposition. 
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1. Introduction 

Distributed quantum computing [I] and quantum cryptographic protocols [2] require 
the transmission of quantum information, or entanglement, over large distances. The 
natural candidate for a transmission line of this kind is optical fibre, with qubits 
encoded in the states of broadband single-photon wavepackets. The ability to transfer 
such 'flying' qubits to a material system, where they could be stored or manipulated 
in a controlled way, would greatly facilitate the implementation of many powerful 
protocols f° r quantum information processing. 

Here we present a scheme for a quantum memory, in which an ultrashort single- 
photon signal pulse is coherently absorbed within an atomic vapour through a two- 
photon Raman transition (see figure ^) . This idea has been explored previously for 
the case of narrowband photons 0, drawing on the dynamical theory developed by 
Raymer and Mostowski |7j. By analysing the problem in terms of its fundamental 
mode structure, we show that the main difficulty with this previous proposal arose 
from poor modematching, because of the use of narrowband fields. More recently, this 
theory was revisited 00 in an extension of the entanglement generation protocol 
proposed by Duan, Cirac, Lukin and Zoller ^Old- ^ n this paper we investigate a 
related but distinct dynamics, which affords the ability to store, rather than generate, 
an entangled photon. In this sense our scheme is similar to that of Fleischhauer et 
al. EH i based on Electromagnetically Induced Transparency. However, in our 

scheme there is no group velocity reduction for the signal photon. This implies that the 
scheme is best suited to short photon wavepackets whose length in space is less than 
that of the atomic medium. The storage does not rely on the resonant phenomenon of 
quantum interference, and this allows for the mapping of temporally short, broadband 
wavepackets into the memory. In our proposal the spatio-temporal structure of the 
signal photon is transferred by a strong classical control pulse to a long-lived collective 
atomic excitation |12| , which we call a spin wave. This is consistent with the fact that 
the internal atomic degree of freedom is often, but not always, an angular momentum. 
After some delay, the stored photon may be read out by sending in another control 
pulse; the frequency and structure of the recovered photon depends upon the shape 
of this readout pulse. The control pulse shapes which optimise the memory efficiency 
may be generalised for resonant quantum memory schemes |14| . 




control pulse 




Figure 1. Left: the level structure of the atoms comprising a quantum memory 
for broadband photons, with bandwidth <5. Right: a schematic of the readin 
process for the quantum memory. 
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The results in this paper are one dimensional, in the sense that they apply to an 
ensemble which is much longer, along the propagation direction, than it is wide (with 
a Fresnel number |15j of order unity or less). A fully three dimensional model, capable 
of investigating the memory evolution with non-colinear signal and control fields, as 
well as the limits imposed on memory fidelity by spontaneous emission and motional 
decoherence, is the subject of a parallel publication JHj- In the one-dimensional limit, 
we show that the full propagation problem is simplified greatly by the use of an 
appropriate modal decomposition. 

In [8] it is shown that the Raman interaction of a CW driving field with an atomic 
vapour can be decomposed as a set of two-mode squeezing transformations between 
light and matter. These modes characterise both the temporal structure of the Stokes 
light produced, and the spatial structure of the correlated spin wave generated in 
the ensemble. Recent work [H] has verified the exact nature of this squeezed mode 
decomposition, shedding light on the properties of such multimode interactions. In 
particular the authors show how canonical evolution of the field operators fixes the 
relationship between the input and output modes for the problem. This analysis of 
Stokes generation is closely connected to the description of the quantum memory we 
present here. We show in this paper that the structure of the dynamical equations 
describing both systems guarantees canonical evolution of the optical and spin-wave 
field operators. Their commutation relations are automatically preserved by the 
solution and we do not need to impose this property as an extra condition. That 
is, the mode structure is intrinsic to the classical electrodynamics of the problem, 
rather than being intrinsically quantum mechanical. We also show that, if group 
velocity dispersion can be neglected, the evolution of the quantum memory subject to 
a control pulse of arbitrary shape is related to the simple case of a CW control field by 
a coordinate transformation. The evolution has a natural decomposition in terms of 
a set of mutually orthogonal input and output modes, and these modes are found to 
be related by a time-reversal-like symmetry. Each mode is associated with a singular 
value which quantifies the efficiency of the transfer from the input light-field mode 
to the stored spin-wave mode. The magnitudes of the singular values depend only 
upon a single parameter which represents the strength of the coupling of the quantum 
memory to the incident signal photon. 



The control pulse, together with the ensemble of three-level atoms, may be treated as 
a two-level system, whose absorptive properties are determined by the control field. 
In [7] a fully quantum mechanical treatment for the spontaneous initiation of Raman 
scattering in a one-dimensional atomic ensemble is introduced. The arguments may 
be adapted to describe the propagation of a signal photon through the memory cell. 

The i th atom in the ensemble with energy eigenstates \ at position Zi and time 
t, couples to the electromagnetic field via the dipole interaction Hamiltonian (setting 



2. Model 



K= 1): 




(1) 



The transition projection operators are defined by 



o-) k {t = Q)=\j i ){k i \, 



(2) 
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in the Heisenberg picture. The electric field E(t, z) — E s (t, z) + E c (t, z) is made up of 
a classical control field with central frequency io c and polarisation unit vector e c , 

E c (t, z) = e c E c £(t, z) exp [-iw c (t - z/c)] + c.c, (3) 

and a single photon signal field with polarisation unit vector e s , written in quantized 
notation as 

/■oc 

E s (t,z) = e s i / g(u))a(t, u>) exp[icoz/c] dco + h.c. (4) 
Jo 

The spatio-temporal dependence of the control field is contained in the dimensionless 
envelope function £(t, z), which is scaled so that its maximum value is 1. E c is then the 
peak amplitude of the control field. In the Heisenberg picture, the mode annihilation 
operators satisfy the usual equal-time commutation relations 

[a(t, Lu),a f (t, lu')] = 5{uj - J). (5) 



The mode amplitude is given by g(u>) = i7T ^ gAc , where A is the cross-sectional area 

of the region within the ensemble excited by the optical fields, c is the speed of light 
and eo is the permittivity of free space J7| . We have also defined the transition dipolc 
moments by 

A) k ee (JM^h), (6) 

where d l is the electric dipole moment operator for the i th atom. 

The signal pulse amplitude is represented by the 'slowly varying' operator 

1 f°° 

A(t, z) = exp [iiv s (t — z/v s )\ / a(t, lo) exp [iwz/w s ] dui, (7) 



/2vr 

where lo s is the centre frequency of the pulse, and v s is the group velocity of the signal 
field. The spin wave to which the signal excitation is transferred is represented by the 
corresponding operator |SJ 

B (t, z ) = z -f= Y] 0"i3 {t) exp I -i (u) s - u c )t- ( — - — ) z & \ 
ozynAp L L \ v s v cJ }) 

with v c the group velocity of the control field. Here the index (3 runs over all the 
atoms lying within a thin slice of the ensemble of thickness Sz, centred at position 
z. n is the atomic number density. The operators A and B act as continuous- mode 
bosonic annihilation operators for the signal pulse and atomic ensemble respectively; 
they satisfy the canonical commutation relations 

[A(t,0),At(t',0)]= S(t-t'), (9) 

[B(0,z),B\0,z')}=S(z-z'), (10) 

where the Dirac delta function 8(x) is the derivative of the Heaviside step function 

, , f 1 for x > 0; 
for x < 0. 



,(8) 
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3. Propagation 

Provided that the common detuning A of the signal and control pulses from single 
photon resonance is much larger than the signal photon bandwidth 5, the excited 
state \rrii) can be adiabatically eliminated. If the ensemble is prepared in the collective 
groundstate |0) = i and if the population of the metastable state |3j) is assumed 
to remain negligible, a linear theory can be used. The Maxwell-Bloch equations, in 
the slowly varying envelope approximation, are then found to be: 

d t B(t,z) = -j*£*(t,z)A(t,z), (12) 

and 



d z + -d t 



A(t,z) = jS(t,z)B(t,z). (13) 
The coupling constant 7 is given, to first order, by 



T .ygs^, ^ (14) 

with d s , d c the dominant transition dipolc moments: 

d s = d l ml ■ e s ; d c = d l m3 ■ e c , Vi. (15) 

The dephasing of the material excitation B can be modelled by appending a decay 
term and a Langevin operator JHj to l|12fl . but here for simplicity we neglect the effects 
of decoherence|: this remains a good approximation if the transverse coherence decay 
time of the spin wave is much larger than the duration of the signal pulse. We define 
the local time r = t — z/v s , whence we obtain 

d T B(T,z) = -7*£*(t- kz)A(t,z)] (16) 

d z A(r,z)= jE(t-kz)B(t,z). (17) 

Here we have expressed the functional form of the control pulse in terms of the 
dispersivity k, given by 

k= — -— . (18) 

V c V s 

In the case that dispersion can be neglected, so that v c sa v s , the control field becomes 
a function of r only, and a Laplace transform over the spatial variable z yields an 
analytic solution to the coupled equations (|16I17I) 7 . This is the situation with which 
the bulk of this paper is concerned. However, even if v c =/= v s , the solutions can be 
obtained numerically In the general case, for arbitrary k, they have the following 
form: 

A(r,L)= f C a (t,t')A(t',0)<1t'+ [ S A (T,z')B(0,z')dz', (19) 
Jo Jo 

B{T,z)= ( C B {z,z')B(0,z')dz' - / S b {z 7 t')A(t', 0)dr', (20) 



Jo Jo 
where L is the length of the ensemble, and T is the duration of the readin process. 
The integral kernels Ca.b, Sa,b are Green's functions which propagate the boundary 
conditions A(t,0) and B(0, z). Because we write the solutions at the output of the 
ensemble, after the pulses have fully traveled through it, these are input-output, or 
scattering, relations. 

% We have also neglected the Stark shift of the levels and |3;}. The former, induced by the 
presence of a single photon, is truly negligible; the latter is larger, but the control pulse centre 
frequency can always be chirped appropriately to cancel its effect. 
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4. Unitarity and Mode Decomposition 

4-1- Quantum memory dynamics 

The coordinates (r, z) may be discretized over an arbitrarily fine mesh, and then the 
above equations (|19I20|) are represented (as in [Sj) to any chosen accuracy by the 
matrix equations 

A L = CaAq + SaBq; (21) 
B T = C B B — S B A a , (22) 

with column vectors Aq,Bq replacing the initial signal and spin- wave amplitudes, 
and with Al,Bt replacing the final signal and spin- wave amplitudes. The integral 
kernels are replaced by the matrices C a,b and Sa,b- This transformation should 
be unitary, since it describes the evolution of quantum mechanical operators, but 
it is not immediately obvious how this is guaranteed by the equations (|16I17|I . As 
with an ordinary lossless beamsplitter, it is conservation of flux which constrains the 
solution. To see this, consider the following continuity relation implied by the evolution 
equations: 

d z A^(r, z)A{t, z) + d T B\r, z)B{t, z) = 0. (23) 
Integrating this expression over a rectangle in (r, z)-space yields the expression 

f T A\¥, z)A(t', z) dr'+ f" B\t, z')B(t, z') dz' = F t (t)+F z {z), (24) 
o Jo 

where F z and F T are functions determined by the boundary conditions. Successively 

setting (r, z) = (0, 0), (0, L), (T, 0), (T, L), we derive the flux-excitation conservation 

condition 

A\ ■ Al + Bj, ■ Bt = Ay ■ Aq + Bq • B Q , (25) 

which must hold for arbitrary initial amplitudes A ,B . Here the 'dot' denotes the 
scalar product of two vectors: A B = A T B, where the superscript T indicates matrix 
transposition. The dagger denotes Hcrmitian conjugation of operators. When applied 
to an ordinary matrix, or a complex number, it is equivalent to complex conjugation. 
In order to avoid confusion, we do not use a dagger to indicate the composition of 
transposition with complex conjugation. We can cast the condition (|25|l in a more 
compact form: 

X f • X = Xl ■ X , (26) 
where we have defined 

X = ( * L ) = UX , (27) 



T 



with 



' B 



H26JI then tells us that the transformation Xq — > X preserves the norm of Xq, and this 
fixes U as unitary. The evolution of the operators describing the memory interaction 
is therefore canonical, as we would expect. This fact has important implications. In 
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particular, multiplying out the relations U T *U = UU T * = I provides us with the 
conditions 

C T /C A + S T B *S B = I; (29) 

C t b *Cb + SI*Sa=I, (30) 

along with a pair of antinormally ordered counterparts. We can now apply the Bloch- 
Messiah reduction ^H] to Q29J1 . For Ij29(l . this consists in spectral decomposition 
|2()| of the positive Hermitian matrix products C t a *Ca and S b *Sb- These must 
commute with one another, and therefore they are both rendered diagonal in the 
same orthonormal basis. Similar conclusions follow for the remaining equations f (|3L)[l 
and antinormal versions). As is shown in [5], the analysis reveals simple relationships 
between the singular value decompositions of the integral kernels, as follows: 



C a (t,t')= Y,H-r)^i{r'), (31) 



S b (t',z)= YM*)Wi(/)> (32) 



i=l 
oo 

E 

i=l 
oo 



C B (z,z') = S2<Pi( z )VtXt(z'), (33) 



S a (t,z')= ^^(t)A 4X4 (z'), (34) 
i=i 

where the functions {4>i\, {V'i}; {x-i}y each form a complete orthonormal basis, 

and where Xi^fii are real, positive singular values for which 

A? +tf = l Vi. (35) 

Equations (|31I35|) . when substituted into l|19l2(Jfl . imply a set of independent beam- 
splitter transformations of ensemble modes and light-field modes. 

4-2. Application to Stokes scattering 

Finally, we note that the above arguments apply in slightly altered form to the analysis 
of stimulated Stokes scattering presented in [§] . In that case, Stokes photons and spin 
wave excitations are generated in pairs, and accordingly it is the flux difference that 
is conserved by the Maxwell-Bloch equations (where now the operator A represents 
the amplitude of the Stokes field): 

d T B(T,z) =7*£*(T-Kz)A t (r,z); (36) 
d z A(r, z) = 7 £(t - kz)B^(t, z). (37) 

The resulting flux-excitation conditions are altered slightly, corresponding to swapping 
the plus sign appearing in (12311 for a minus sign. We can write them in the form 

X n ZX = XI T ZX Q ; (38) 
X T ZX ] = X T ZXl, (39) 

where 
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is the z-projection Pauli spin matrix. The solution to (|36I37|) is now written 

X = CX + SXl (40) 
where the matrices C,S are of the form 

C A \ „ _ I S A 



Cs { o c B )- Ss [-s B »)■ W 

Substitution of this ansatz into the sum of i|38|) and (|39|) provides us with the conditions 
C T *ZC + S T ZS* = I: (42) 
C r *ZS + S T ZC* = 0. (43) 

These relations fix the inverse transformation to (|4UII as 

X = ZC T *X - ZS T X\ (44) 

and substitution of this into the flux conditions (|38I39[) yields the antinormal condition 

CC T * - SS T * = Z. (45) 

Inspection of the components of the equations 142145(1 reveals all the conditions 
required of the matrices Ca,b,Sa,b to apply the Bloch-Messiah reduction. This 
then yields a mode decomposition for the Stokes scattering problem, which is given 
explicitly in j^j and used to identify the dynamics with that of the well-known 
two-mode squeezer. The property of separability which allows the dynamics to be 
understood in terms of independent modes is therefore also related to a form of flux 
conservation. 

We now show that when dispersion can be neglected, the relations (|31H34|) . 
describing the memory, simplify even further. Under this approximation, only one set 
of functions {4>i}, which may be computed by numerical diagonalisation of a known 
function, is involved in the dynamics. This approximation is consistent with the 
requirement of adiabaticity S <C A for a broadband signal photon, since the refractive 
index of the ensemble varies slowly with frequency far from resonance. 

5. Dispersionless case 

We return to the equations i|ltjll7|) . and set k = 0, so that the ^-dependence of the 
control field £ vanishes. We introduce the dimensionless 'control pulse area' 

e = e (r)^ 7 y|^ T |f(r')| 2 dr', (46) 

where the factor E = J Q T \£(t')\ 2 dr' appears in order to normalise this new variable 



so that it takes values between and C = jy/LE. We also introduce a normalised 
spatial variable £ = 72 \J E / L. Again, £ runs from to C. For simplicity we have set 
7 to be real; alternatively any complex phase in 7 can always be incorporated into the 
field £. In this new coordinate system the equations ll(jll7fl become 

d c a = P ; d e (3 = -a, (47) 

where we have defined the dimensionless annihilation operators a(e, C) = 
VEA(t,z)/£(t), and [3{e,() = VlB(t,z). We have now eliminated the variation 
of the control field from the dynamical equations; the problem is reduced to that of 
a CW control field coupled to the ensemble for a brief time. The solutions for the 
signal amplitude ac(e) = a(e, C) at the exit face of the ensemble, and the spin wave 
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amplitude /3c (C) = P(G, C) at the end of the readin process, are given in terms of the 
initial conditions a(e, £ = 0) = ao(e) and /3(e = 0, Q = Po(C) by the expressions 

a c (e)= [ [G 1 (e-x,C)a {x)+G (C-x,e)Po(x)}dx, (48) 



(i 



/3c(C)= / [t?i(C-a:,C0A)(a;)-Go(C-a:,Oao(a;)] ds. (49) 
Jo 

There are now just two integral kernels, given by 

G (p,q)= J (2VW), (50) 

G 1 (p,q) = 6{p)-Q(p)^J 1 (2^pq), (51) 

V p 

with the n th Bessel function of the first kind denoted by J„, and where the Heaviside 
step function ensures that the convolutions in Ij48l49|l respect causality. We write 
the matrix representation of the solution I|48I49|) as 

a c = G 1 a Q + G /3 ; (52) 

p c = Gi(3 - G <* . (53) 

Note that both of the matrices G\ and Go are per symmetric. That is, they are both 
symmetric about their main anti-diagonal (or skew diagonal). To see this, observe that 
the integral kernel G\ (as it appears in 1)480*) varies only as a function of the quantity 
e — x. That is to say, its contours are all parallel to the line e = x, which corresponds 
to the main diagonal of G±. G\ is therefore invariant under reflection about its main 
anti-diagonal. The corresponding symmetry for Go follows from the Hcrmiticity of 
the kernel Go(p,q) in its arguments. These symmetries allow us to decompose the 
kernels using input and output modes related by time reversal (or equivalently space 
reversal) : 



G (C -x,e)=J2 ^(e)Ai^(C - x), (54) 

oo 

Gi(e -x,C)=J2 <f>i(t)»i<f>i(C ~ x), (55) 

i=l 

where there is now just a single set of real modefunctions {&} such that 

cC 

4>i(x)4>j (x) Ax = 5ij . (56) 
6. State transfer 

Following the treatment in |Hj we define output mode operators using the mode 
expansions 

a c (e)= (57) 

i 

c (C)=J2^OBi. (58) 
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For the input operators, we use the time-reversed expansions 

a Q {e)=Y J MC-e)a u (59) 

i 

f3 (0= J2<t>i(C-0k- (60) 

i 

The solutions I|48I49I) can now be written in the form 

Ai = /^a,, + Xibi , (61) 
Bi = Hib{ - Xiai , (62) 

with 

[At, A}} = [Bi, B)] = [a u a]] = [h, b}] = 5 ly (63) 
Equations (|61I62I) show that the interaction of a signal photon with the memory 
cell can be viewed as a beamsplitter transformation on a mode-by-mode basis. 
Maximising the fidelity of the state transfer then amounts to minimising the terms 
Hidi. The magnitudes of the fit are determined by the coupling parameter C; a large 
value corresponds to an optically thick ensemble — a high absorption. In figure |21 
the singular values of the matrix Go are plotted as a function of C. These are found 
by multiplying (|54"|) through by the i th input mode and integrating, using We 
obtain the following relation: 
rC 

G (C - y, x)0 4 (C -y)dy = X^x), (64) 



corresponding to the simple eigenvalue equation 

J (2y/xy)<l>i(y) dy = Xi<f>i(x), (65) 

which we solve numerically using a 500 by 500 square grid. We see that the interaction 
is dominated by the lowest mode for small C, but as C increases, higher modes become 
significantly coupled. Modematching is therefore particularly important when the 
interaction is weak, as is the case for large detuning. Near-complete state transfer for 
the lowest mode (Ai ~ 1) can be achieved for an ensemble with a coupling strength 
C> 2. 



7. Modematching 

The efficiency of the memory is greatest for the lowest mode (j>\ . Therefore we require 
that the shape of this mode in real time matches the temporal structure of the signal 
photon wavepacket we send into our memory cell. The relationship between the shapes 
of the normalised real time input modes $j(r), and the mode shapes as a function of 
control pulse area e, is determined by the temporal shape of the control field. The 
connection is straightforward, 

Mr)^^£(r)MC~e(r)}, (66) 

but not analytically invertible, so a simple optimisation was performed in order to find 
the control pulse shape which matches a given input photon to the memory mode. 
The result for a signal photon with a Gaussian input profile, matched to the mode 
4>\(C — e) with the highest coupling, and the simplest shape, is shown in figure for 
C = 2. The photon is almost completely absorbed, although a small portion remains, 
due to the limitations of the modematching optimisation. 
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1 2 3 4 5 

c 

Figure 2. The five largest singular values of the kernel Go, plotted as a function 
of the coupling parameter C. The singular values are necessarily positive; for 
clarity they are plotted both above and below the origin. 




Figure 3. Top: the intensity (A*(r. z)A(r. z)) of a Gaussian signal photon as 
it propagates through an atomic ensemble with C = 2. That is, a Fock state 
= a^|vac), where = J^° x £(t)A(t, 0) dr destroys a photon with a wave 
packet amplitude £(r) oc exp {-21n2[(r - t )/ct] 2 }, with a = 1/8, r = T/2 and 
T = L = 1. |vac) denotes the vacuum. The absorption is nearly complete; it is 
limited by the quality of the modematching optimisation. Bottom: the optimised 
control field intensity (blue/black) is shown alongside the initial signal field 
intensity (green/grey). The latter has been scaled for clarity. The lowest mode 
(corresponding to the largest singular value Ai) was used in the optimisation. 



8. Readout 

Once a properly modematched photon has been read in to the quantum memory, the 
ensemble is left in the output mode </>i(£), with probability amplitude Ai. We now 
consider the effect of sending a second control pulse, propagating in the same direction 
as the initial control pulse, into the ensemble. The centre frequency, bandwidth and 
intensity of this readout pulse may differ from that of the first control pulse (herein 
the readin pulse). We can use all of the results developed thus far to analyse the 
interaction of the readout pulse with the ensemble, but we should be careful to keep 
track of any parameters which differ between the readin and readout stages. Let 
us use a superscript r to indicate those quantities associated with the readout. If 
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there is no significant change in the control and signal group velocities at readout 
(so that v T c w v r s sa v c w u s ), then the readout spin wave annihilation operator B r 
is phasematched to the readin spin wave operator B. That is, since the Stokes shift 
uj s — uj c — uj r s — uj r c is fixed, the phase factor appearing in (JSJ) is unchanged at readout, 
and we have B r (0, z) = B(T, z). Here again we have neglected any decoherence of the 
spin wave over the storage period. This provides us with one boundary condition; the 
second is that the signal field begins in its vacuum state at the start of the readout 
process: 

(A rt (r,0)A r (r,0)) =0, Vr. (67) 

At the end of the readout process, some portion of the stored excitation has been 
transferred back to the optical field. The efficiency of the readout depends upon the 
degree to which the spin wave mode 4>\{Cz/ L) (written in terms of the ordinary spatial 
variable z) overlaps with the input modes for the readout process, which are of the 
form $1 [C r (l — z/L)]. The functions {4%} solve the readout eigenvalue equation 

j (2V5fl#(y)dy = (68) 

Each readout input mode </>^ \C r {\ — z/L)] is transferred to the optical field with 
amplitude A£", according to the relation (|61(l . with these amplitudes set by the size 
of the readout coupling parameter C r . A measure of the fidelity of the memory is 
the expectation value of the output photon number operator Af, which is just the 
probability of retrieving a photon from the ensemble at readout, given that a single 
modematched photon was sent in with the readin pulse. This retrieval probability is 
given by 

Af= [ {A r \r 1 L)A r (t, L)) dr 



= -L f (a r \e r ,C r )a r (e r ,C r ))de r 

i—l 
i=l 



i 



fh (69) 



where in the penultimate step we have used the transformation l|tjl|) and the boundary 
condition l|H7|) (setting (a^a-) = (&[ a£) = 0), and where we have defined the i th 
overlap 

fi [ L ft [C r (l - z/L)} MCz/L) dz. (70) 



L 

In figure 01 the variation of Af is plotted as a function of the readin and readout 
coupling parameters C and C' r . We observe that the retrieval probability is low if 
either C or C r is small, as we expect from a weakly driven interaction. However, even 
if the photon is almost completely absorbed at readin, as is the case for the region 
C > 2, the probability of recovering the stored photon approaches unity only slowly as 
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C r is increased. Furthermore, as the readin process is driven harder (as C rises above 
2), the retrieval probability falls. These observations are explained by considering the 
overlaps /, of the readout modes with the spin wave mode. For small C, <pi(Cz/L) 
is monotonic, and relatively flat. It therefore has a large overlap /i with the lowest 
readout mode, which for moderate values of C r is essentially the mirror image of the 
spin wave mode. As C increases, the spin wave mode becomes more asymmetric, so 
that its overlap with the lowest readout mode falls. It is then necessary to pump 
the readout process harder, so that higher readout modes, with which the spin wave 
overlaps significantly, are efficiently coupled to the optical field. The readout efficiency 
therefore grows as the readin coupling C is decreased. But as C becomes small, so also 
does the probability that the photon was stored during the readin process at all. The 
retrieval probability is maximised along the line C ~ 2, which represents the optimal 
coupling for the readin process. However, a readout coupling parameter in excess of 
10 is required to achieve AT > 0.95. 

The time reversal symmetry between the input and output modes makes the 
readout for this scheme a non-trivial problem. Simply repeating the readin process 
(so that C — C) results in poor performance of the memory. However, the dramatic 
increase in coupling strength required to extract the stored excitation fully, may make 
a naive increase in control pulse energy at readout prohibitively difficult to realise. An 
alternative method to boost the coupling is by significantly reducing the bandwidth 
of the readout pulse, along with its detuning A r . The photon recovered from the 
memory in this way would be frequency-shifted (according to the Raman resonance 
condition) , and temporally stretched (since its bandwidth would be greatly diminished 
as well). The temporal profile of the emitted photon would then be determined by 
the shape of the readout pulse, as well as the shapes of the output modes involved 
in the readout process. A 'photon transducer' with the ability to store broadband 
photons and convert them into narrowband photons with controllable frequency and 
shape could prove to be a useful tool for optical quantum information processing. 

Finally, we comment that the symmetry of the modes suggests the possibility 
of reading the stored excitation using a readout pulse propagating in the reverse 
direction. Switching the propagation direction of the readout pulse would send 
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4>i(Cz/L) — > fa [C(l — z/L)], so that the spin wave mode would overlap exactly with 
the lowest readout mode with C = C r . Unfortunately, a counter propagating readout 
pulse would not have sufficient momentum to convert the spin wave into a counter 
propagating anti-Stokes photon. Another way to say this is that the readin and 
readout spin wave operators are not phasematched in this situation; the operator 
B r acquires a rapidly spatially varying phase factor, so that the overlap integrals fi 
vanish. However, an implementation of this scheme in the solid state (for example, 
replacing the atomic ensemble with an ensemble of semiconductor charge quantum 
dots |2U), might allow for the possibility of controlling the ensemble level structure 
with external fields. Then this kind of reverse readout might be achieved with the 
use of guasz-phasematching 22 , in which the sign of the readout coupling parameter 
C r is periodically flipped along the length of the ensemble. If the wavevector of this 
modulation matches the size of the wavevector mismatch between the readin and 
readout spin waves, the readout photon can be retrieved. Alternatively, an adiabatic 
reversal of the sign of the Stokes shift over the whole ensemble, between the readin 
and readout stages (so that the energies of the states |3i) and |li) are swapped at 
readout), would enable the phasematching of a high frequency counter propagating 
readout pulse, with the spin wave transferred to a counter propagating photon with 
lower frequency. Exploration of the possibilities for coherent manipulation of collective 
quanta in gaseous and solid-state ensembles will form the basis of continuing research. 



9. Transverse profile 

Here we briefly show that the modes {</>,;} also describe the transverse profile 
of the signal field. If the operators a and f3 are allowed to depend upon the 
radial displacement r from the z-axis, the Maxwell-Bloch equations, in the paraxial 
approximation 123f , take the modified form 

^K + id C a(p,C,e)=i/3(p,C,e); (71) 

d e P(p,C,e) =-a(p,C,e), (72) 

where p is a normalised position vector for a point in the plane orthogonal to the 
propagation direction, a radial distance r = yj2Lv s /uj s p from the z-axis. The extra 
term Vp is a 2-dimensional Laplacian in the dimensionless cylindrical polar coordinates 
parameterising the point p. Here we have retained the assumption that the control 
field £ has no radial dependence over the area A. The dimensionless radial variable 
p takes values from to p max = ^/io s j (2Lv s )r max , where 7rr max = A. The Fresnel 
number T is given by |15| 



2nv s J L 

— Pmax- (73) 

The limit of the paraxial approximation is usually defined by T = 1, so we set 

Pmax 1 • 

The solution to the equations H71I72|) is found by alternately Fourier transforming 
the variables (p, C) an d (P; e ) E3| 123 • We find: 

a c (p,e)= [ d 2 p' [ dx[R(p-p',C)G 1 (e-x,C)a (p',x) + 

J A JO 



Pc(p,C) 
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R(p -ff,C- x)G (C - x, e)p (p', x)] ; (74) 
f d 2 p' f dx[R(p-p'X-x)G 1 (C-x,C)(3 (p',x)- 

JA JO 

R(p - p\ C)Go(C - x, C)ao(p', a:)] , (75) 
where d 2 p' denotes an integral over a circular patch with unit radius, and where 
the diffraction kernel is given, up to a normalisation constant N, by 

R{p,x) = Nexp(-iC\p\ 2 /x)/x. (76) 

This kernel is singular for vanishing x, since the paraxial approximation breaks down 
in the near-field. However, in the limit that the interaction is negligible close to the 
exit face of the ensemble, we can make the replacement R(p,x) — > R(p,C). The 
diffraction kernel then factorises out of the e and C integrals in I|74I75|) . We introduce 
the spectral decomposition 

oo 

R(p-p',C) =X>(P)<^(P'), (77) 

3=0 

where the o~j are real eigenvalues and {(fj} is a complete orthonormal set of paraxial 
modes, satisfying 

/ fj(pM{p)d i p = 6 jk . (78) 

J A 

We assume that the spatial profile of the signal field depends only upon the radial 
coordinate p. Cylindrical symmetry of the boundary conditions then allows us to 
neglect modes with any dependence on azimuthal angle. These cylindrically symmetric 
modes are then given by 

cp j (p)=e i » 2 <f> 1 j (p% (79) 

where the superscript 1 indicates that the function </>j solves the eigenvalue equation 
(|65|l with C set equal to 1. We define input and output mode operators, as in 
H57I58I59I60|I . by the relations: 

OO 

ao(p,e)= ^(C-^GoHi; ( 80 ) 

(j.i)=o 

OO 

/3o( P ,C)ee MC - Ofj(j>)ba , (81) 



and 



"c(p,e)= J2 M £ )<Pj(p)Aii 5 ( 82 ) 

(i,j)=G 

OO 

Pc{p,0= MC)<PjO>)Bij , (83) 



with 



[^ii> A U = [Ba,B)d\ = 4il = N. fo Ll = ( 84 ) 
We then obtain the beamsplitter transformations: 

Aij = o-j (piOij + Xibi 3 ) ; (85) 

Bij = o-j {Hihj - AiCiij) . (86) 
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The lowest mode is dominantly coupled, with o~i — 0.995. The shape of this dominant 
spatial mode is well approximated by a Gaussian with beam waist w s = 1.45 (in 
the normalised units of the variable p). The control should be slightly less tightly 
focussed, with a waist w c > 3w s , to justify the approximation that the control field 
may be treated as a train of plane waves. If such appropriately focussed Gaussian 
beams are used for the signal and control fields, only the lowest paraxial mode is 
involved in the interaction. The index j may be dropped from the relations 1)85186(1 , 
and setting o~\ ~ 1, we recover the transformations 1)61162(1 . and the one-dimensional 
treatment is justified. 

10. Experiment 

We are currently developing an implementation of this scheme in atomic thallium 
vapour. The gross structure of thallium has a A-type triplet, where the excited ls 2 S 1 /2 
state is coupled strongly to the ground state &p 2 P\/2 and the metastable state 6p 2 P 3 / 2 
(with an electric-dipole- forbidden radiative lifetime of 0.16 s). The splitting of these 
last two levels is very large (around 8000 cm -1 ), so that the signal and control pulses 
remain spectrally distinct, even when their bandwidths are large. The more stringent 
limit on the bandwidth of the photon that can be stored comes from the adiabaticity 
condition: in order that no real population is transferred to the excited state, atomic 
collisions must occur at a negligible rate and the signal bandwidth 6 should be 
much smaller than the detuning A. The interaction cross-section scales as the inverse 
square of this detuning, so a compromise must be made between the bandwidth S that 
is stored and the length L of the ensemble. We hope to demonstrate the storage and 
retrieval of ultraviolet pulses (~ 390 nm) with S w 10 nm and pulse lengths of around 
100 fs. 

11. Summary 

We have shown that it is possible to map a broadband photon wavepacket into 
an atomic spin wave using a Raman interaction. The relationship between the 
input state, the storage state and the output is clarified by the introduction of a 
mode decomposition. The modes are expressed in terms of a single universal set of 
eigenfunctions of a particular integral kernel. The form of these eigenfunctions depends 
only on the coupling strength of the memory interaction, and not upon the shapes of 
the pulses. 
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